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We discuss a balanced homodyne detection scheme with 
imperfect detectors in the framework of the operational ap- 
proach to quantum measurement. We show that a realistic 
homodyne measurement is described by a family of opera- 
tional observables that depends on the experimental setup, 
rather than a single field quadrature operator. We find an 
explicit form of this family, which fully characterizes the ex- 
perimental device and is independent of a specific state of the 
measured system. We also derive operational homodyne ob- 
servables for the setup with a random phase, which has been 
recently applied in an ultrafast measurement of the photon 
statistics of a pulsed diode laser. The operational formula- 
tion directly gives the relation between the detected noise 
and the intrinsic quantum fluctuations of the measured field. 
We demonstrate this on two examples: the operational uncer- 
tainty relation for the field quadratures, and the homodyne 
detection of suppressed fluctuations in photon statistics. 

PACS Number(s): 42.50.Dv, 03.65.Bz 



I. INTRODUCTION 

Homodyne detection is a well known technique in de- 
tecting phase-dependent properties of optical radiation. 
In quantum optics it has been widely used in studies and 
applications of squeezed light |ffl] . A statistical distribu- 
tion of the outcomes of a homodyne detector has recently 
found novel applications in the measurement of the quan- 
tum state of light via optical homodyne tomography ||] 
and the direct probing of quantum phase space by pho- 
ton counting || . The phase-sensitivity of homodyne de- 
tection is achieved by performing a superposition of the 
signal field with a coherent local oscillator by means of 
a beamsplitter It was an important observation jjj 
that in a balanced scheme, with a 50%:50% beamsplitter, 
the local oscillator noise can be cancelled by subtracting 
the photocurrents of the detectors facing two outgoing 
beams. Then, in the limit of a classical local oscillator, 
the statistics of difference photocounts is simply a sig- 
nal quadrature distribution rescaled by the amplitude of 
the local oscillator Hf7|. Therefore, balanced homodyne 
detection is an optical realization of an abstract quan- 
tum mechanical measurement of the field quadratures 
described by a quantum observable xg. The statistical 



outcomes of an ideal measurement of xg\xg) = xg\xg) 
are described by the spectral measure 



p(xg) = ( \xg)(xg\ [ 



(1) 



Although the spectral measure contains all the relevant 
statistical information about the homodyne measure- 
ment, it corresponds to a quantity that is measured by 
an ideal noise-free detector. Due to this property xg will 
be called an intrinsic homodyne quantum observable. 

However analysis of the homodyne setup with imper- 
fect detectors JtJ shows that the relation between the 
statistics of the difference counts and the quadrature 
spectral distribution is in fact more complicated. The 
distribution measured in a real experiment is smoothed 
by a convolution with a Gaussian function of width de- 
pendent on the detector efficiency. Consequently, realis- 
tic homodyne detection cannot be straightforwardly in- 
terpreted as a measurement of the intrinsic field quadra- 
tures Xg. 

A recent experimental application of homodyne detec- 
tion to the reconstruction of the photon number distri- 
bution of a weak field from a pulsed diode laser Q has 
shown that a homodyne setup with the fluctuating phase 
9 is a powerful tool in measuring phase-insensitive prop- 
erties of light. However, it is not possible to associate 
with this setup any spectral measure even in the case of 
perfect detectors. Therefore, homodyne detection with 
the random phase cannot be described in terms of mea- 
suring any intrinsic quantum observable. 

It is the purpose of this paper to show that homodyne 
detection provides an interesting and nontrivial exam- 
ple of a realistic quantum measurement leading to op- 
erational quantum observables, i.e., to quantum opera- 
tors that depend on properties of a specific experimental 
setup used in the homodyne detection. In particular, 
these operational observables will depend on the detec- 
tor losses described by a quantum efficiency n and on the 
phase 9 of the local oscillator used to probe the signal 
field. Such operational observables provide a natural link 
between the quantum formalism and raw data recorded 
in a realistic homodyne experiment. 

General features of the operational approach, with ref- 
erences to earlier literature are given in IJ. The main 
conclusions of this approach, if applied to the homodyne 
measurement, can be summarized as follows. A quantity 
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delivered by the home-dyne experiment is a propensity 
density Pr(a) of a certain classical variable a. This den- 
sity is given by an expectation value of an a-dependent 
positive operator valued measure (POVM), denoted by 
H(a): 



Pr(a) = (W(a)). 



(2) 



Thus, the POVM given by H (a) corresponds to a realistic 
homodyne detection and is the mathematical representa- 
tion of the device dependent measurement. In one way of 
looking at quantum measurements, the emphasis is put 
on the construction and properties of such POVMs. In 
such an approach, in realistic homodyne detection, the 
spectral decomposition dxg\xg)(xg\ of the intrinsic ob- 
servable ig, is effectively replaced by the POVM daft (a). 
Consequently the moments of Pr(a) can be represented 



as 



= J da o n Pr(o) 



H 



(3) 



defining in this way a family of operational homodyne 
quantum observables 



n 



daa n H(a) 



(4) 



where the index TL stands for the homodyne detection 
scheme associated with the given POVM. This family 
characterizes the experimental device and is independent 
on a specific state of the measured system. 

In this paper we derive and discuss the family of op- 
erational observables for balanced homodyne detection 
with imperfect photodetectors. We show that for bal- 
anced homodyne detection an exact reconstruction of 
the POVM 7i(a) and of the corresponding operational 

quantum quadratures be performed. Thus, 

homodyne detection provides a nontrivial measurement 
scheme for which an exact derivation of the correspond- 
ing POVM and the operational observables is possible. 
The interest in construction of this operational algebra 
is due to the fact that the number of physical examples 
where the operational description can be found explicitly 
is very limited jlC|j . We show that the algebraic prop- 
erties of the Xg\j differ significantly from those of the 

powers of xg. In particular x g 2 ^ n ^ (i^ ' H ) 2 - This prop- 
erty will have immediate consequences in the discussion 
of the uncertainty relation with imperfect detectors. 

This paper has the following structure. First, in Sec. ||, 
we derive the POVM and the generating operator for the 
operational observables. Their explicit form is found in 
the limit of a classical local oscillator in Sec. III. Given 



this result, we discuss the operational uncertainty rela- 
tion in Sec. [V. In Sec. [v| we derive the family of opera- 
tional homodyne observables for the homodyne detector 
with a random phase between the signal and the local 
oscillator fields, and relate them to the intrinsic photon 



number operator. These calculations link the homodyne 
noise with fluctuations of the photon statistics, and can 
be useful in the time-resolved measurement of the prop- 
erties of pulsed diode lasers. Finally, Sec. [Vl] summarizes 
the results. 



II. GENERATING OPERATOR FOR 
HOMODYNE DETECTION 

The family of the operational homodyne quantum ob- 
servables defined in Eq. (Q) can be written conveniently 
with the help of the generating operator 



dae lXa H{a). 



(5) 



Operational quantum observables are given by deriva- 
tives of the generating operator at A = 0: 



,{n) 



1 d" 

i n d\ n 



(6) 



A=0 



This compact representation will noticeably simplify fur- 
ther calculations. 

We will start the calculations by finding the generat- 
ing operator for the homodyne detector. In a balanced 
setup, the signal field described by an annihilation opera- 
tor a, is superimposed on a local oscillator b by means of 
a 50%:50% beamsplitter. The annihilation operators of 
the outgoing modes are given, up to the irrelevant phase 
factors, by the relation 



1 

V2 



1 1 

1 -1 



(7) 



We will assume that the local oscillator is in a coherent 
state |/3) lo- If another state of the local oscillator is con- 
sidered, our formulae can be generalized in a straightfor- 
ward manner by averaging the results over an appropriate 
Glauber's P-representation. 

A quantity recorded in the experiment is the statistics 
of the difference counts between photodetectors facing 
the modes c and d. The difference of the counts AN 
corresponds to the classical variable, denoted before as 
a, recorded in a homodyne detection experiment. The 
POVM 'H(AN) describing this detection scheme can be 
easily derived. It is clear that this POVM is an operator 
acting in the Hilbert space of the signal mode. Its explicit 
form can be found with the help of standard theory of 
photodetection pJ] : 

H(AN) = ^lo{\P)(P\lo 

ni-?i2-A]V 
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where r\ is the quantum efficiency, assumed to be iden- 
tical for both the detectors. In this formula the partial 
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trace is over the local oscillator mode and a marginal av- 
erage with a fixed value of AN is performed. We will 
now convert this POVM into the generating operator ac- 
cording to Eq. (|5|). As we discuss later, the POVM H 
and consequently the generating operator Z-j-c have their 
natural parametrization, independent on the LO inten- 
sity. Before we find this scaling, we will use £ instead of 
A as a parameter of the generating operator: 

oo 

AjV=— oo 

= Tr LO {|/3)(/3| LO 

: exp^e* - l)c f c + rj(e _i4 - l)Sd] : }. (9) 

Let us transform this expression to the form which 
does not contain the normal ordering symbol. For this 
purpose we will use the technique developed by Yuen 
and Shapiro Q consisting of extending the Hilbert space 
by two additional modes c„ and d v and constructing the 
fields annihilation operators 



c<j = \fri'c + V 1 VCv, d d = ^/rjd + y/l - r]d v . (10) 

The generating operator can be written in the extended 
four-mode space using these operators as 



ZH(£) = T±LO,v{\0){P\l>Ot 

: exp[(e l « - l)c\c d + (e^ - l)4<i d ] : }, (11) 

where Tr v denotes the trace over both the vacuum modes 
c v and d v . We can now apply the relation jl2[ 



exp[(e^ — l)w^«] : = exp(i£v'v) 



(12) 



valid for an arbitrary bosonic annihilation operator v, 
which finally gives: 

= Tr LOi!J {|/3)(/3| LO <g> |0)(0|„ exp^(4c d - d\d d )]}. 

(13) 

This expression contains the most compact form of the 
homodyne POVM. The exponent in Eq. ( |l~ij ) resembles 
the one from Eq. (^), with c, d replaced by Cd, dd and 
the detector efficiency equal to one. It is known that 
there is a physical picture behind this similarity. An im- 
perfect photodetector can be equivalently described by 
an ideal detector preceded by a beamsplitter with the 
power transmissivity equal to the quantum efficiency of 
the real detector, assuming that the vacuum state enters 
through the unused port of the beamsplitter. Mathe- 
matically, this construction corresponds to the so called 
Naimark extension of the POVM into a projective mea- 
sure on a larger Hilbert space |l(J . 



III. APPROXIMATION OF A CLASSICAL 
LOCAL OSCILLATOR 

When the local oscillator is in a strong coherent state, 
the bosonic operators b\ b' can be replaced by c-numbers 
(3,(3*. However this approximation violates the bosonic 
commutation relations for the pairs c, c' and d, (V , which 
have been used implicitly several times in the manipu- 
lations involving Z-n(£). Therefore some care should be 
taken when considering the classical limit of the local 
oscillator. 

We will perform the approximation on the exponent of 
Eq. (JH^). We will replace the quantum average over the 
state |/3) by inserting (3, (3* in place of b, w and keep only 
the terms linear in (3. This gives 



c d c d 



V2 



(14) 



The operator in the brackets has a form analogous to 
Eq. ([To]) with the combination (c v +d v )/ \/2 as a vacuum 
mode. Consequently, the imperfectness of the photode- 
tectors in the balanced homodyne detection can be mod- 
elled by superposing the signal on a fictitious vacuum 
mode before superposing it with the local oscillator and 
attenuating the amplitude of the local oscillator field by 
y/rj. This observation has been originally made by Leon- 
hardt and Paul |l3| , and is an example of the Naimark 
extension involving a nonquantized local oscillator. 

Under the approximation of a classical local oscilla- 
tor, it is now easy to perform the trace over the vacuum 
modes with the help of the Baker-Campbell-Hausdorff 
formula. This yields 

Z n (0 = expKVl - V )\(3\ 2 /2} exp^^at + (3*a)}. 

(15) 

The exponent exp[— £ 2 r)(l — r))\{3\ 2 /2] introduces a spe- 
cific ordering of the creation and annihilation operators 
in the generating operator. Therefore the detector ef- 
ficiency r] can be related to the ordering of the oper- 
ational observables. For example, for r\ = 1/2 we get 

Zh{£) = exp(i£/3*a/2) exp(i£/?a 1 '/2), i.e. the generating 
operator is ordered antinormally. 

The expansion of the generating operator into a power 
series of i£ gives 



1 d n Z n 




l-ry 



x 



(16) 



where H n denotes the nth Hermite polynomial and x$ is 
the standard quadrature operator 
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V2 



(17) 



expressed in terms of the creation and annihilation oper- 
ators of the signal field and dependent on the local oscil- 
lator phase 9 defined as [3 = \(3\e lS . In the terminology of 
the operational approach to quantum measurement, xg 
is called an intrinsic quantum observable, since it refers 
to internal properties of the system independent of the 
measuring device || . 

It will be convenient to change the parameter of the 
generating operator in order to make the derivatives ( |l6| ) 
independent of the amplitude of the local oscillator. A 
scaling factor which can be directly obtained from an 
experiment is the square root of the intensity of the local 
oscillator field measured by the photodetector. We will 
multiply it by \[2 in order to get the intrinsic quadrature 
oper ator (p"7j ) in the limit rj — ► 1. Thus substituting £ = 
A/ y/2r]\(3\ 2 yields the generating operator independent of 
the amplitude of the local oscillator: 



IV. OPERATIONAL UNCERTAINTY RELATION 

With explicit forms of operational observables in hand, 
we can now analyze their relation to the intrinsic quadra- 
ture operator. For this purpose, let us look at the first 
lowest-order operational quadrature observables: 



x 



(!) _„l/2 



H 



;.(2) 



x e'n- r J[ x e+ 2r? 

i(3) _ 



1 — rj 



•3/9 / 31 — 7] 

x e -' n = r\i [x g + -—xe 



2 V 



(21) 



Z n {\, 9) = exp[-A 2 (l - j])/4 + iXy/^/2(, 



>*;>t 



e-*»a)]. 



The imperfectness of photodetectors influences the oper- 
ational observables in two ways. The first one is a trivial 
rescaling of the observables by the powers of y/rj, the sec- 
ond way is a contribution of the lower-order terms to 
the operational counterparts of Xg. In order to see its 
conseque nces let us investigate the rescaled operational 
variance (AA) 2 - (AN) 2 : 



(18) 



Sxl n 



1 



{AN) : 



(22) 



The derivatives of Zu{\ 
ily of the operational observables £ 
detector: 



ive the final form of the fam 

(n) 



n 



V 2i 



H,, 



From the definitions of the operational operators it is 
n for the homodyne clear that this operational variance involves xp u and 

x^ n . The combination of these two operators is in gen- 
eral different from the intrinsic variance. Because of this 
the operational dispersion of xe is: 



(19) 



The algebraic properties of the operational observables 
are quite complicated, since ij' is not simply an nth 

power of Xg 1 ^ n . Thus a single operator does not suf- 
fice to describe the homodyne detection with imperfect 
detectors. Complete characterization of the setup is pro- 
vided by the whole family of operational observables. In 
fact the operators x 6 , w define an infinite algebra of op- 
erational homodyne observables for an arbitrary state of 
the signal mode. As mentioned above, for rj = 50%, the 
general formula reduces to antinormally ordered powers 
of the intrinsic quadrature operators: 



Sx 2 en 



x(2) 



w 



H' 



Axi 



2q 



(23) 



where Axe — y/ (stg) — (xe) 2 is the intrinsic quantum dis- 
persion of the quadrature xe ■ This intrinsic dispersion is 
enhanced by a term coming from the imperfectness of the 
detectors. Thus, the imperfectness of the photodetectors 
introduces an additional noise to the measurement and 
deteriorates its resolution. 

Using the above result we can derive the operational 
uncertainty relation for the quadratures related to the 
angles 9 and 8' 



1 



H 



2™/2 



\{xbY 



(20) 



5xe n & x 6'H > V ( Aa^Acc 



1 — 7] 

2r] 



(24) 



This expression shows that the operational operators are 
in some sense equivalent to a prescription of ordering of 
the intrinsic operators. This prescription is dynamical in 
character, i.e., it depends on the efficiency rj of the detec- 
tors used in the homodyne detection. In fact the homo- 
dyne operational algebra is defined by a one-parameter 
family of dynamical orderings defined by the generating 
operator derived in this section. 



Again, an additional term is added to the intrinsic un- 
certainty product. This situation is similar to that 
in Rcf. |l4j where it was argued that taking into ac- 
count the measuring device raises the minimum limit 
for the uncertainty product. However, that discussion 
concerned a simultaneous measurement of canonically 
conjugate variables, which is not the case in homo- 
dyne detection. Using the intrinsic uncertainty relation 
AxgAxe> > | sm(8—8')\/2 we get the result that the right 
hand side in the operational relation ( p4] ) is not smaller 
than (?7| sin(0 - 8')\ + 1 - n)/2. 
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One may wonder if the definition of squeezing is af- 
fected by the operational operators. Let us consider the 
two quadratures 8xen an d $ x e+f H - In this case the op- 
erational uncertainty, 



8xg H Sxe+ 



2H 



1 

> -. 



(25) 



is independent of 77. However, it has to be kept in mind 
that only a part of the operational dispersion comes from 
the field fluctuations. The easiest way to discuss this is 
to rewrite Eq. (|23|) to the form 



1 
71 



(26) 



which shows that the operational dispersion is a 
quadratic average of the intrinsic field dispersion Ax$ 
and the detector noise l/v2 that corresponds to the vac- 
uum fluctuation level. These contributions enter with the 
weights 77 and 1 — ry, respectively. Therefore if a squeezed 
quadrature is measured with imperfect detectors, the ob- 
served dispersion is larger than the intrinsic one. 



V. HOMODYNE DETECTION WITH RANDOM 
PHASE 

Homodyne detection is used primarily to detect phase- 
dependent properties of light. However, it has been re- 
cently shown that even a setup with a random phase 
between the signal and local oscillator fields can be a 
useful tool in optical experiments ||. Although in this 
case the phase sensitivity is lost, the homodyne detector 
can be applied to measure phase-independent quantities 
and such a setup presents some advantages over a sin- 
gle photodetector. First, the information on the statis- 
tics of the field is carried by the photocurrent difference 
between the two rather intense fields. Within existing 
detector technology, this quantity can be measured with 
a significantly better efficiency than the weak field itself. 
Secondly, the spatio-temporal mode that is actually mea- 
sured by the homodyne detector is defined by the shape 
of the local oscillator field. Consequently, application of 
the local oscillator in the form of a short pulse allows the 
measurement to be performed with an ultrafast sampling 
time. This technique has been used in Ref. § to measure 
the time resolved photon number statistics from a diode 
laser operating below threshold. The achieved sampling 
time was significantly shorter than those of previously 
used methods. 

The photon number distribution and other phase- 
independent quantities are reconstructed from the aver- 
age of the random phase homodyne statistics calculated 
with the so-called pattern functions |l5| ] . For commonly 
used quantities, such as the diagonal elements of the den- 
sity matrix in the Fock basis, these pattern functions 
take a quite complicated form. In this section we will 



consider observables that are related to the experimental 
data in the most direct way, the moments of the homo- 
dyne statistics with randomized phase. We will derive 
the family of operational observables and relate them to 
the powers of the photon number operator n — a)a. 

The generating operator for homodyne detection with 
random phase Z-r. {7Z stands for the random phase) is 
obtained readily from Z-h by averaging it over the phase 
6. This gives 



ZkW 



d0 
2^ 



Z n (\,6) = e- A2 / 4 : Jo (\^/2r,a^aj 



(27) 



where Jo is the Bessel function of the Oth order. With the 
help of the result derived in the Appendix, the normally 
ordered form of the Bessel function can be transformed 
into the following expression: 



Ztz(X) = e-* 2 / 4 L A (r,\ 2 /2), 



(28) 



where the index of the Laguerre polynomial is the the 
photon number operator. The Laguerre polynomials 
with a operator valued index is defined by the decom- 
position in the Fock basis. 

The family of operational observables is given by the 
derivatives of the generating operator 



'n 



1 d n 
^dA" 



(29) 



A=0 



Since the homodyne statistics averaged over the phase is 
even, the odd derivatives disappear. A straightforward 
calculation yields the operators for even 11 = 2m: 



(2m - 1)!! 



E 



(2rj) k 



k\ 



k=0 

xfiifi — 1) . . . (n — k + 1). 



(30) 



This formula shows that xj^ is a polynomial of n 
of the order of m. Therefore the first m moments of 
the photon number distribution can be computed from 



~(2) 



(2m) 



The two lowest order observables are 



given explicitly by 



,(2) 



r\n 



a(4) _ 
TZ ~~ 



- \r 1 2 h 2 +ri(2-i 1 )n+- 



(31) 



It is seen that even in the case of ideal noise- free detectors 



v ti 7^ i^n') 2 an d the family of the operational observ- 
ables has nontrivial algebraic properties. Inversion of the 
above equations yields: 
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x iz 2 



(32) 



As an illustration, let us express the normalized photon 
number variance Q = ((n 2 )-(n) 2 -(n))/(n) @ in terms 

of the expectation values of x^) and x$ . This variance 
is used to characterize the sub-Poissonian statistics of 
light. After some simple algebra we arrive at 



1 5 {Xfo ) 



V 



(33) 



Thus, the variance Q can be read out from the two lowest 
moments of the homodyne statistics with the randomized 
phase. The photodetector efficiency r\ enters into the 
above formula only as an overall scaling factor. This 
result is analogous to that obtained for the setup with a 
single imperfect detector, and is due to the fact that Q 
describes normally ordered field fluctuations. 



VI. CONCLUSIONS 

We have presented the operational description of the 
balanced homodyne detection scheme with imperfect 
photodetectors. For homodyne detection it is possible 
to derive exact expressions for the POVM and the cor- 
responding algebra of operational operators. The result 
of these calculations shows that a whole family of oper- 
ational observables rather than a single operator should 
to be used to discuss a realistic setup. This family allows 
one to easily relate the experimentally observed fluctua- 
tions to the intrinsic properties of the system. 



ACKNOWLEDGMENTS 



ta : =]T(f + £ )»<n| 

n=0 

n=0 k=0 V 7 
oo oo / \ 



A;— n—k 



If we now assume the convention that 



(Al) 



for 



n < k, the range of the second sum can be extended to 
k = to oo. ft is then natural to denote it as a binomial 
coefficient o f the operator n. Comparing the equal powers 
of e in Eq. (Al) yields a very compact representation of 
the normally ordered powers of aJa in terms of h: 



{a)) k a k = M 



(A2) 



Expanding the normally ordered Bessel function in 
Eq. (E7T) and applying the above identity gives: 



: J (Xy^2r]¥a\ : = ^ 

fc=0 

oo 

= £ 



k=0 



(fc!) 2 (a ) a 
k\ 



i#i(r?A 2 /2), 



(A3) 



where the Laguerre polynomial with the photon number 
operator index 



Ln(x) = 22 L n (x)\n)(n\ 



(A4) 



is defined analogously to the binomial coefficient via de- 
composition in the Fock basis. 
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APPENDIX: 

In this Appendix we present details of the transfor- 
mation of the generating operator Z-jz(X) into the form 
that does not contain the normal ordering symbol. Let 
us start by rewriting Eq. (12) to the form exp(ea ji a) = 
(1 + e) n and decomposing its right hand side of in the 
Fock basis {|n)}: 
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